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sorbed.  In  contrast  to  previous  calculations,  it  is  also  found  that  the  ordinary  mode  is  well  absorbed 
at  oblique  incidence  to  the  toroidal  magnetic  field.  This  finding  indicates  the  feasibility  of  heating 
devices  of  higher  densities  than  would  otherwise  be  possible. 
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THEORY  OF  SECOND  HARMONIC  ELECTRON  CYCLOTRON 
RESONANCE  HEATING  OF  TOKAMAK  PLASMA 

I.  INTRODUCTION 

It  has  been  shown  that  electron  cyclotron  resonance  heating  at  the  fundamental  harmonic 
could  be  very  promising  for  heating  plasma  to  ignition  for  future  tokamak  reactors.1-4  In  this 
paper,  we  shall  examine  the  characteristics  of  electron  cyclotron  resonance  heating  at  <o  —  2 
w«(0)  using  a  fully  relativistic  ray  tracing  code  in  toroidal  geometry.  Here  o>  is  the  wave  fre¬ 
quency  and  to„(0)  is  the  electron  cyclotron  frequency  at  the  center  of  the  plasma.  If  the 
plasma  is  in  thermal  equilibrium,  the  results  presented  here  could  also  be  applied  to  the  study 
of  plasma  emission  at  o>  —  2o>„(0)  by  invoking  KirkolTs  Law. 

The  motivation  for  considering  second  harmonic  heating  is  that  the  accessibility  condi¬ 
tions  allow  the  heating  of  higher  density  tokamaks  with  the  second  harmonic.  For  example,  for 
wave  vector  angles  not  too  far  from  normal  incidence  to  the  magnetic  field,  the  accessibility 
conditions  at  w  —  <u„(0)  are  o>f2,(0)  >  o>£,(0)  for  the  ordinary  mode  denotes  the  electron 
plasma  frequency),  and  2«2,(0)  >  w^(0)  for  the  extraordinary  mode  launched  from  the  high 
magnetic  field  side  of  the  tokamak  (waves  launched  from  the  low  magnetic  field  side  are  inac¬ 
cessible).  On  the  other  hand,  for  heating  at  w  —  2««(0),  the  above  accessibility  conditions  (for 
near  normal  incidence)  become  4w2,(0)  >  <o£(0)  for  the  ordinary  wave  and  2a><?,(0)  >  to^(0) 
for  the  extraordinary  wave  (see  Sec.  II).  Note  that  in  the  extraordinary  wave  case  the  largest 
allowable  density  is  the  same  for  <o  -  w„(0)  and  &>  —  2tuc,(0),  but  that  the  latter  case  has  the 
advantage  that  waves  can  be  launched  from  either  the  high  or  the  low  field  sides  of  the 
tokamak  plasma  column.  More  significantly,  for  the  ordinary  wave  at  u  -  2to„(0)  the  largest 
allowable  desnity  is  four  times  that  for  the  ordinary  wave  at  o>  -  to„(0)  and  twice  that  for 
extraordinary  waves.  To  put  this  situation  in  perspective  we  can  express  the  plasma  0(0)  in  the 
center  of  the  plasma  in  terms  of  a>^(0)/a>c2f(0) 
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0(0)  -  (8%)  (r/10KcV)  {«^(0)/ft>f,(0)], 

where  T  is  the  plasma  temperature  (assumed  equal  for  electrons  and  ions).  For  example,  at 
T  “  lOKeV,  the  maximum  allowable  0(0)  is  8%,  16%,  and  32%  for  the  ordinary  wave  at 
to  -  «„(0),  the  extraordinary  wave,  and  the  ordinary  wave  at  to  -  2o>«(0),  respectively.  (Note 
that  the  plasma  0  averaged  over  the  tokamak  cross-section  would  typically  be  0(O)/3].  The 
possibility  of  heating  high  density  (i.e.  high  0)  tokamaks  evidently  depends  on  whether  the 


ordinary  wave  at  to  ■■  2«„(0)  can  be  effective  absorbed  in  the  plasma  interior.  A  major  result 
of  the  present  study  is  that  this  is  indeed  found  to  be  the  case  (Secs.  Ill  and  IV).  This  result  is 


based  on  the  consideration  of  oblique  angles  of  incidence  to  the  magnetic  field  and  contrasts 


with  previous  results  obtained  for  normal  incidence 
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II.  ACCESSIBILITY 

For  simplicity,  in  this  section  we  consider  a  straight  plasma  cylinder  with  density 
N(r)  ( dN/dr  <  0  and  N(a)  s  0),  and  magnetic  field  B  —  B(r,9)  z0.  To  simulate  the  toroidal 
vacuum  field  dependence  on  the  mtuor  radius  we  take  B  -  B0  (1  -  x/R0)~ '  where 

• -  f  t 

x  —  —  r  cos  0  and  R0  is  the  major  radius  of  the  torus.  Furthermore,  we  assume  incident  waves 
launched  at  9  —  0  or  6  —  ir  with  no  9  component  of  the  incident  wavenumber  ( k9  *■  0).  Since 
B  has  no  9  component,  kt  remains  zero  as  the  wave  propagates  inward.  (Note  that  these 
simplifications  are  not  made  in  our  ray  tracing  calculations  to  be  presented  in  Sec.  III.) 

To  find  the  accessibility  condition  at  to  —  2cd[(,(0),  we  proceed  as  in  Ref.  1  and  examine 
the  cold  plasma  dispersion  relation.  We  obtain  three  possible  reflection  points, 


*1-1 

(1) 

X2  -  (1  -  S2)  (1  +  Y) 

(2) 

X3  -  (1  -  S2)  (1  -  Y) 

(3) 

X\  is  an  ordinary  wave  reflection  point,  X}  is  an  extraordinary  wave  reflection  point,  and  X2  is 
an  ordinary  [extraordinary]  wave  reflection  point  for  S2  >  Y(Y  +  2)-1  [for  S2  <  Y  (Y  +  2)-1]. 
Resonance  for  the  extraordinary  mode  occurs  at 

X,  s  1  -  Y1  (4) 

and  mode  conversion  between  the  ordinary  and  extraordinary  modes  occurs  at 

Xc  =  1  +  Y2  -  ~  f  -)2 
452 

where  X  =  otjf/u2,  Y  —  o>ce/<o,  S  -  sin©  «•  ck||/u>,  0  is  the  angle  of  incidence  (for  normal 
incidence  5-0)  and  k||  denote  the  component  of  k  parallel  to  B. 

In  order  for  the  wave  to  be  accessible  to  the  center  of  the  plasma,  there  must  be  no 
reflection,  resonance,  or  mode  conversion  points  between  the  plasma  edge  and  the  center  of  the 
plasma  (where  maximum  absorption  occurs).  By  examination  of  the  expressions  for  Xi,  X2, 
XJt  Xr,  and  Xc,  we  find  that  this  can  be  achieved  for  F(0)  -  1/2  if 

*(0)  ^  (1  -  S2)/ 2,  (6) 
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for  the  extraordinary  mode,  and  if 


J  1  forS2  ^  1/3 

Jf(0)  *  1  3/2(1  -  S2)  for  S 2  ?  1/3  (7) 

for  the  ordinary  mode.  These  accessibility  conditions  are  displayed  graphically  in  Figs.  1  and  2 


Note  that,  in  contrast  to  heating  at  the  fundamental  harmonic,  the  extraordinary  mode  is  acces 


(0),  which  is 


For  the  extraordinary  mode,  the  highest  density  accessible  is  ml 


no  improvement  over  heating  at  w  -  w„(0).  For  the  ordinary  mode,  the  density  limit  is 


(0),  which  is  probably  sufficient  for  any  future  tokamak 
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III.  NUMERICAL  CALCULATION  OF  LINEAR  ABSORPTION 

To  find  out  how  efficient  the  heating  at  the  second  harmonic  is,  we  use  a  relativistic  ray 
tracing  code  previously  developed  and  described  in  Ref.  1  to  study  the  linear  absorption  of  a 
single  pass  of  both  the  ordinary  and  the  extraordinary  modes  in  toroidal  geometry.  The  density 
and  temperature  profiles  are  assumed  to  be  parabolic. 

Two  types  of  tokamaks  are  studied:  (A)  "reactor  size",  an  examples  is  the  UWMAK  III 
(minor  radius  —  270  cm.,  mqjor  radius  —  810  cm);  and  (B)  "present-day  size",  an  example  is 
the  Princeton  Large  Torus  (minor  radius  —  40  cm.,  major  radius  -  132  cm.).  The  results  of 
the  calculations  are  summarized  in  Table  I  and  II. 

We  see  that  the  extraordinary  mode  is  always  efficiently  absorbed  for  the  range  of  angles 
considered,  i.e.,  close  to  perpendicular  incidence.  This  high  absorption  agrees  with  previous 
theoretical  predictions.2  3 

More  importantly,  our  results  show  that  the  second  harmonic  ordinary  mode  is  also  quite 
efficiently  absorbed  (especially  for  "reactor-size  tokamaks")  at  oblique  incidence  (an  estimate  of 
the  ordinary  mode  absorption  for  slab  geometry  at  oblique  incidence  is  presented  in  the  next 
section).  Even  at  exactly  normal  incidence,  the  second  harmonic  ordinary  mode  absorption  is 
larger  than  previous  predicted2-3  (probably  due  to  a  combination  of  ray  wandering  and  relativis¬ 
tic  effects  which  are  not  included  in  the  analyses  of  Refs.  2,3). 
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IV.  CALCULATION  OF  THE  ORDINARY  WAVE  TRANSMISSION  COEFFICIENT 

The  absorption  of  wave  energy  at  the  second  harmonic  occurs  mainly  as  a  result  of  the 
coupling  of  the  electron  gyromotion  with  the  circular  component  of  the  wave  field  that  rotates 
in  the  electron  direction.  For  this  reason  the  damping  of  a  perpendicularly  propagating  ordinary 
wave,  for  which  the  perpendicular  electric  field  is  zero,  is  much  smaller  than  that  of  the 
extraordinary  mode.  However,  if  the  ordinary  wave  is  propagating  at  some  finite  angle  with 
respect  to  the  perpendicular,  it  acquires  a  perpendicular  component  of  the  electric  field;  and 
therefore,  it  becomes  more  highly  damped. 

To  calculate  this  damping  we  evaluate  the  transmission  coefficient  for  propagation  through 
ansi  -  2o>„  resonant  layer. 


T  =  exp  (-  2irt)) 

where  ttt)  -  J"  dx  k,(x)  is  the  optical  depth, 

k,(x)  -  Im  (E0+  .  D,(A:0,x)  ■  E0)  (E0+  •  D 0(kB)  •  (8) 

aK0 

is  the  imaginary  part  of  the  local  wave  vector2,  E0  and  k0  are  the  wave  electric  field  and  wave 
vector  as  determined  from  the  cold  plasma  dielectric  tensor  D0,  and  E>i(/c0,  x)  is  the  imaginary 
part  of  the  correction  to  D„  resulting  from  the  finite  Larmor  radius  effect  at  the  second  har¬ 
monic. 


To  determine  k,(x)  we  assume  k-B  <<  |fc||£|  so  that  the  perpendicular  component  of 
the  electric  field  is  small  compared  with  the  parallel  component.  We  then  find2 

i-i 


E*  •  Di  •  E„  -  -rl£_l2— ?■ 


1  I  r-  15 


«2  mt( O2 


x 

y  -I-  /  0 


and, 


E/  •  D0  •  E0  *•  y|£2 1  k0c2/<o}e 

where  £||  is  the  parallel  component  of  the  wave  electric  field,  £_  is  the  relevant  circularly 
polarized  component  of  the  perpendicular  electric  field,  and  L  is  the  scale  length  of  the  mag- 
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inhomogeneity,  L  1  =  \  VB/B\.  The  quantity  £_  can  be  calculated  in  terms  of  E„  from  the 
cold  plasma  dielectric  tensor6 


1  kuk.  k 

i£_/£,.|2ai  *4. 


,  k2c2 

-K,  KrK,  -  i-V  (Kr  +  K,) 


where  K,  -  1  -  w£/(6«£),  K,  -  1  - 


and  k2c2/at 2  -  1  -  o>£/(4a>2f).  Clearly  as  — •  0,  £_  — »  0.  Upon  performing  the  integration 

of  *,(x)  over  x  we  find. 


-I*0|£_|2/|£M|2. 


-2  2  01  1  '  l  ii l  * 

2o>ZemecL 

which  applies  for  kjk0  <1.  In  general  we  must  add  to  Eq.  (9)  the  small  contribution  that  oc¬ 
curs  for  exactly  perpendicular  propagation2.  However,  the  contribution  given  in  Eq.  (9)  will  be 
larger  than  this  term  by  a  factor  of  order  (kjk0)2  \  Tjmec2]~x. 


A  slab  model  numerical  code  was  used  to  test  Eq.  (9).  The  analytical  and  numerical 
results  were  found  to  be  in  excellent  agreement  (<  2%  discrepancy)  with  8  <  20°,  i.e.  when 
k n  «  k0.  One  notes  from  Eq.  (9),  however,  that  7j(a/f|f)  is  a  very  sensitive  function  of  the 
angle  8.  In  a  tokamak  because  of  the  complicated  geometrical  factors,  this  angle  varies  consid¬ 
erably  from  point  to  point  in  the  path  of  the  wave.  Consequently,  while  Eq.  (9)  serves  to 
demonstrate  the  possibility  of  strong  ordinary  wave  absorption  at  the  second  cyclotron  har¬ 
monic,  it  has  limited  applicability  to  systems  where  kH  is  not  conserved.  This  indicates  the 
necessity  of  utilizing  the  ray  tracing  code  as  we  have  reported  in  the  previous  section. 
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V.  CONCLUSION  AND  DISCUSSION 

Generally  speaking,  if  the  plasma  temperature  is  above  1  keV,  both  the  ordinary  and  the 
extraordinary  modes  can  be  efficiently  absorbed  at  either  the  first  or  the  second  cyclotron  har¬ 
monic.  For  Te  >  IkeV,  the  determining  factor  in  regard  to  the  choice  of  mode  is  the  ratio  of 
(0^(0)  to  <uc.e(0).  Table  III  summarizes  the  possible  methods  of  heating  in  various  regimes  of 
ca^fOi/to^fO).  In  the  regime  wpe(0)  <  a>((,(0),  ordinary  mode  heating  at  the  first  cyclotron 
harmonic  frequency  is  the  most  desirable  because  it  allows  wave  injection  from  the  low  field 
side.  In  the  regime  cu^fO)  <  ai^fO)  <  V2caIf(0),  ordinary  mode  can  no  longer  be  employed 
and  there  is  a  trade-off  between  the  extraordinary  mode  at  the  first  harmonic  and  the  extraordi¬ 
nary  or  ordinary  mode  at  the  second  harmonic.  Either  choice  has  a  disadvantage,  namely,  the 
former  mode  has  to  be  launched  from  the  high  field  side  while  the  latter  modes  require  a 
higher  wave  frequency.  In  the  regime  V2eoC(.(0)  <  oj^CO)  <  2a>„,(0),  the  only  choice  is  the 
ordinary  mode  at  the  second  harmonic.  From  the  above  consideration,  it  is  clear  that  electron 
cyclotron  resonance  heating  at  the  second  harmonic  frequency  is  a  viable  and  possibly  the  only 
alternative  for  reactor  size  tokamaks  which  generally  have  iope(0)  >  <d„,(0).  The  main 
difficulty  involved  in  second  cyclotron  harmonic  heating  is  the  development  of  a  high  frequency 
wave  source.  In  the  case  of  UWMAK  III,  the  required  wave  frequency  is  225  GHz.  In  view  of 
recent  progress  made  in  the  gyrotron  research,  however,  this  does  not  appear  to  be  an  insur¬ 
mountable  problem.  For  example,  a  100  GHz,  1  MW  gyrotron  with  30%  efficiency  has  already 
been  reported7  and  an  experiment  aimed  at  240  GHz,  20  kW,  and  25%  efficiency  is  currently 
under  way  at  the  Naval  Research  Laboratory.  Thus,  there  is  a  reasonable  prospect  for  obtaining 
the  required  wave  source  in  the  next  generation  of  gyrotrons. 

Relatively  large  emission  of  the  second  harmonic  ordinary  and  extraordinary  mode  have 
been  observed.8  9  The  emission  level  of  the  second  harmonic  ordinary  mode  is  much  larger 
than  would  be  predicted  on  the  basis  of  previous  theoretical  results,  which  only  consider  normal 
incidence  and  neglect  relativistic  effects.  On  the  other  hand,  our  results  are  consistent  with 
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Fig.  2  —  Maximum  accessible  value  of  ot^/ to3,  in  the  center  of  the 
plasma  versus  S3  for  the  extraordinary  mode. 
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Table  I  —  Transmission  Characteristics  of  "Reactor- 
size"  tokamak.  The  magnetic  field  and  electron  tem¬ 
perature  in  the  center  are  B0  —  4.05  T.  and  T0  -  2keV. 
N0  denotes  the  density  in  the  center  of  the  device. 


Mode 

0 

(0) /to  ce,( 0) 

Tran. 

Coeff. 

N0(cm  3) 

mm 

0° 

1.9 

■7YTHB 

5.77  x  10l4 

mm 

10° 

1.9 

5.77  x  1014 

02 

20° 

1.9 

5.77  x  1014 

02 

0° 

1.3 

0.53 

2.7  x  1014 

02 

10° 

1.3 

0.0 

2.7  x  1014 

02 

20° 

1.3 

0.0 

2.7  x  1014 

X2 

10° 

1.3 

0.0 

2.7  x  1014 

X2 

20° 

1.3 

0.0 

2.7  x  1014 

Table  II  —  Transmission  Characteristics 
of  "Present-day  size"  tokamak. 
^(0)/o*„(0)  -  0.71,  T„  -  2keV, 
B„  -  3.2  T.,  and  N0  -  5  x  1013cm  ~3. 


Mode 

0 

Tran 

Coeff. 

02 

10° 

0.68 

02 

20° 

0.60 

X2 

5° 

0.0 

X2 

30° 

0.0 

Table  III  —  Methods  for  electron  cyclotron  resonance  heating  in 
various  plasma  density  regimes 


Regime 

Mode/Harmonic 

Low  Field  Side 

High  Field  Side 

Injection 

Injection 

to„(0)  „ 

to„(0) 

01 

yes 

yes 

perpendicular 

XI 

no 

oblique 

02 

yes 

oblique 

X2 

yes 

perpendicular 

A 

e  e 

5  * 

o  © 

A 

XI 

no 

yes 

oblique 

02 

yes 

yes 

oblique 

X2 

yes 

yes 

perpendicular 

rM 

V 

V  . — . 

©  © 

3  3 

V 

02 

yes 

yes 

oblique 

